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Abstract 



Final-state interactions in the response of a many-body system to an exter- 
nal probe delivering large momentum are normally described using the eikonal 
approximation, for the trajectory of the struck particle, and the frozen ap- 
^ ' proximation, for the positions of the spectators. We propose a generalization 

of this scheme, in which the initial momentum of the struck particle is explic- 
itly taken into account. Numerical calculations of the nuclear matter response 
at 1 < |q| < 2 GeV/c show that the inclusion of this momentum dependence 
leads to a sizable effect in the low energy tail. Possible implications for the 

analysis of existing electron-nucleus scattering data are discussed. 
PACS numbers: 13.6.Le, 25.30Fj, 25.30Rw 
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I. INTRODUCTION 



Final state interactions (FSI) of fast nucleons produced in electron-nucleus scattering 
at large momentum transfer have long been known to exert a significant effect on the co- 
incidence (e, e'p) cross section; moreover, they provide most of the strength observed in 
the low-energy loss tail of the inclusive (e, e') cross section (see, e.g., Ref. [|T|). The main 
effect of FSI is a damping of the motion of the struck particle, which can be qualitatively 
described in terms of the imaginary part of the nuclear optical potential. However, since 
nucleons in nuclei are strongly correlated, one must improve upon the optical potential ap- 
proach rooted in a simple mean-field description of nuclear dynamics, if one hopes to develop 
a fully quantitative treatment of FSI. It is very important to realize that nucleon-nucleon 
(NN) correlations, leading to large density fluctuations and to the appearance of high mo- 
mentum components in the nuclear wave function, strongly affect both the initial and the 
final state in electron- nucleus scattering and must be consistently taken into account. 

A theoretical description of (e, e') processes including correlation effects has been de- 
veloped in Ref. (hereafter referred to as I) and successfully employed to analyze (e, e') 
data for momentum transfer in the range 1 GeV/c < |q| < 2 GeV/c 0-^. A similar ap- 
proach, formulated in analogy with the theoretical treatment of FSI in neutron scattering 
from quantum liquids, has been proposed in Ref. [^. 

The treatment of FSI discussed in I, commonly referred to as Correlated Glauber Approx- 
imation (CGA), rests on the assumptions that (i) the struck nucleon moves along a straight 
line with constant velocity (eikonal approximation), and (ii) the spectator nucleons are seen 
by the fast struck particle as a collection of fixed scattering centers (frozen approximation). 
The resulting inclusive cross section can be written as a convolution integral, involving the 
cross section evaluated within the Plane- Wave Impulse Approximation (PWIA), i.e., evalu- 
ated in the absence of FSI, and a folding function embodying FSI effects. The CGA entails 
the same set of approximations as the Glauber theory of high-energy proton scattering off 
nuclei 0, which has been successfully applied for over forty years. 
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The results of I show that FSI produce a huge enhancement of the inclusive cross section 
in the region of oj « ujqe, where ljqe is the energy transfer corresponding to elastic 
scattering off an isolated stationary nucleon. While this enhancement brings theory and 
experiment into agreement over a broad range in u, the calculated cross section substantially 
overestimates the data in the extreme low-cu tail (roughly corresponding to values of the 
Bjorken scaling variable x = {Q"^ /2muj) > 2, where = |qp — uj"^ and m denotes the 
nucleon mass). In order to reproduce the tail of the measured cross sections, the imaginary 
part of the free-space A^A^ scattering amplitude, which determines the shape of the CGA 
folding function, must be modified in such a way as to reduce the effect of FSI. 

As pointed out in I, A^A^ scattering in the nuclear medium may in principle differ 
markedly from scattering in free space. For example, Pauli blocking and dispersive cor- 
rections are known to be important at moderate energies 0. However, their effects on the 
calculated cross sections have been found to be small in the kinematical region spanned 
by the data analyzed in Refs. 0-^]. Corrections to the A^A^ amplitude associated with the 
extrapolation to off-shell energies are also expected to be small 

A different type of modification of the A^A^ cross section, originating from the internal 
structure of the nucleon, may play a more significant role. It has been suggested P,pil| that 



elastic scattering on a nucleon at high momentum transfer can only occur if the nucleon is 
found in the Fock state having the lowest number of constituents, so that the momentum 
transfer can be most effectively shared among them. Within this picture, a nucleon is in a 
very compact configuration after absorbing a large momentum q. It then travels through 
nuclear matter experiencing very little FSI, until its standard size is recovered on a charac- 
teristic time scale that increases with |q|. In the limit of infinite momentum transfer, FSI 
are totally suppressed, and the nuclear medium is said to exhibit color transparency. 

The results of the calculations of Refs. show that inclusion of the effects of color 

transparency according to the model of Ref. |Tl[], with no adjustable parameters, greatly 
improves the agreement between theory and data, yielding a satisfactory description of the 
low-energy loss tail of the nuclear inclusive cross sections for > 1.5 (GeV/c)^. 



To firmly establish the occurrence of color transparency in (e, e') processes, the accu- 
racy of the approximations underlying the CGA must be carefully investigated and either 
validated or transcended. In this paper, we introduce a treatment of FSI which improves 
upon the CGA, in that it allows one to take account of the initial momentum of the struck 
nucleon. Within this approach, the response can no longer be written as a simple convo- 
lution integral. However, it can still be expressed in terms of the spectral function and a 
generalized folding function, in a form displaying explicit dependence on both the initial and 
final momenta of the struck particle. 

The theoretical description of nuclear-matter response is discussed in Section II, where 
we outline the development of a systematic scheme that improves upon the PWIA and 
includes FSI effects. The details of the many-body calculation of the generalized folding 
function within the Fermi-HyperNetted Chain (FHNC) approach are traced in Section III. 
The ensuing numerical results are presented and analyzed in Section IV, with particular 
attention to the generalized folding functions and the nuclear matter y-scaling functions at 
different values of |q|. The final section summarizes our findings and states our conclusions. 

II. NUCLEAR-MATTER RESPONSE 
A. Plane- Wave Impulse Approximation 

The analysis carried out in Refs. |@-|^ required a full calculation of the nuclear cross sec- 
tion, including the electromagnetic vertex, as well as the use of spectral functions adapted 
to finite targets (such as those obtained within a local-density approximation [^]). In ad- 
dition, since the typical momentum transfers lie in the 1-2 GeV/c range, consistent use of 
relativistic kinematics was essential. 

In this paper, we will avoid these complications and focus on the nonrelativistic response 
of infinite nuclear matter to a scalar probe, defined by 
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= y £ e^^-+^o)*(0|pte-^^V.|0) . (1) 

Here H and |0) denote the nuclear Hamiltonian and the corresponding ground state sat- 
isfying the Schrodinger equation H\Q) = Eq\0). The time-dependent density fluctuation 
operator Pq(t) is constructed as 

k 

where aj^ and are nucleon creation and annihilation operators, respectively. Note that 
the definition given in Eq. (|l|) can be readily generalized to describe the electromagnetic 
response by replacing pq with the appropriate current operator. 

Retaining only the incoherent contribution to the response, which is known to be domi- 
nant at large |q|, Eq. (P can be rewritten in the form 

5(q,^) = /^e^(-+^o)*^(q,t), (3) 

with 

S{q,t) = J dRdR'%{R')e''''<{R'\e-'^'\R)e''''''^oiR) , (4) 

where {R} = {ri, r2, r^} specifies the spatial configuration of the A-nucleon system, 
\E'o(-R) = (-R|0) is its ground-state wave function, and the propagator Ua{R, R'',t) = 
(i?'|e~*^*|i?) represents the amplitude for the system to evolve from configuration R to 
configuration R' during the time t. The wave function \l'o(-R) can in principle be evaluated 
within nonrelativistic nuclear many-body theory. On the other hand, the nonrelativistic 
approach cannot be used to obtain the A-particle propagator Ua{R, R'',t), since - in the 
kinematical regime under consideration - the struck nucleon typically carries a momen- 
tum larger than the nucleon mass. In view of the fact that a fully realistic and consistent 
calculation of Ua{R, R'',t) remains intractable, one must resort to simplifying assumptions. 

A systematic approximation scheme can be developed by first decomposing the Hamil- 
tonian according to 
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H = Ha-i + T, + Hi, (5) 

where Hx_i is the nonrelativistic Hamiltonian of the fully interacting [A — l)-particle spec- 
tator system and Ti denotes the Hamiltonian describing a free nucleon. The term 

^^/ = E^i., (6) 

where Vij is the NN potential, accounts for the interactions between the struck particle and 
the spectators. 

The PWIA amounts to setting if/ = in Eq. (1), thus disregarding FSI altogether. The 
resulting A-particle propagator factorizes into the product of the interacting {A — l)-particle 
propagator and the free-space one-body propagator describing the struck nucleon: 



f/i 



WiA(i?, R'; t) = Ua^i{R, R'; t)t/o(ri, r[; t) , (7) 

where {R} = {r2, ...,r^} specifies the configuration of the spectator system. Eq .(0) clearly 
shows that within the PWIA, nuclear dynamics only appears through Ua-i, while the treat- 
ment of the relativistic motion of the struck nucleon reduces to a trivial kinematic problem. 

We next express the PWIA to S'(q, t) (and thereby the response) in terms of the nucleon 
spectral function P(k, E'), which by definition gives the probability of removing a nucleon 
with momentum k from the nuclear ground state, leaving the residual system with excitation 
energy E. Introducing spectral representations for both Ua-i and Uq (see, e.g., Ref. [0), 
we obtain 

^PWiA(q, t) = J-0^sJdEP{p- q, E)e-(^-^°+^'')* , (8) 



which leads to the familiar result 12 



5pwiA(q, ^) = / ^ / dE P(p - q, E)5(^ - E - E,) , (9) 
where Ep = |pp/2m denotes the kinetic energy of a nucleon carrying momentum p. 
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B. Inclusion of Final-State Effects 



In order to improve upon the PWIA, one needs to devise a set of sensible approximations 
to treat the contributions to the A-particle propagator coming from the FSI Hamiltonian 
Hj. As a first step, we make the replacement 

which essentially amounts to assuming that the internal dynamics of the spectator system 
and its FSI with the struck particle do not affect one another, and can therefore be completely 
decoupled. Within this picture, the spectator system evolves during the time t as if there 
were no struck particle moving around, while the fast struck particle sees the spectator 
system as if it were frozen at time t = 0. 

Under this assumption, which implies that the configuration of the spectator system does 
not change due to interactions with the fast struck nucleon, we can use completeness of the 
{A — l)-particle position eigenstates to rewrite the propagator in the simple factorized form 

UA{R,R';t) = J ci^"(^'|e-^^^-i*|^")(r;,^"|e-^(^i+^^)*|ri,^)5='(^-i)(^ - R") 
= Ua^i{R, R'; t) {r[, ^|e~^(^^+^^)*|ri, R) 

= Ua-i{R, R'; t) U,{r,R, v\R- 1) . (11) 

Evaluation of Ui{riR,r[R]t) in general requires a functional integration over the set of 
trajectories ri(r) satisfying the boundary conditions ri(0) = ri and ri(t) = r[ (see, e.g., Ref. 
||13|). However, for large momenta of the struck nucleon, the evaluation can be drastically 
simplified by invoking the eikonal approximation, i.e. by assuming that the particle moves 
along a straight trajectory with constant velocity v = (r'^ — ri)/t, so that ri(r) = ri + vr . 
Within this approximation, the propagator Ui{riR,r[R;t) takes the factorized form [|l^ 

f/i(ri^, v[R; t) = t/o(ri, r'^; t)f/p(ri, R; t) , (12) 

where the eikonal propagator Up{ri, R;t) is given by 
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f/p(ri,i?;t) = exp 

V being the A^A^ potential and f/o(ri, r'^^; t) the free-space nucleon propagator. 

Expanding the exponential appearing in the right-hand side of Eq. (|13|), one obtains a 
series whose terms are associated with processes involving an increasing number of interac- 
tions between the struck nucleon and the spectators. The terms corresponding to repeated 
interactions with the same spectator can be summed up to all orders by replacing the bare 
A^A^ interaction v with the coordinate space t-matrix rq(r), which is related to the A^A^ 
scattering amplitude /q(k) at incident momentum q and momentum transfer k through 

97r r rl^h 

Using the above results together with spectral representations of both f/o(ri, r'^; t) and 
Ua-i{R-, R'', t), the response can finally be expressed as 

5(q, co) = J ^e^(-+^o)t I ^e-^^^* ^ e'^^-^M^Jp - q)Mon(p, q; t) , (15) 

where the sum extends over the {A — l)-particle states satisfying the Schrodinger equations 
HA-i\n) = En\n). We have introduced the definitions 

Mon(k) = J dRe''' '''%{R)<^niR) , (16) 

with $n(-R) = {R\n), and 

Mon(p,q;t) = / rfi?e-^(P-'i) '-^$:(i?)^o(i?)t/p(ri,/?;t) • (17) 

In the limit f/p(ri, R; t) 1, we have Mon(p, q; t) Mon(p — q), and the response S'(q, u) 
given by Eq. (|1^) reduces to the PWIA result. 

It is to be emphasized that the calculation of the response according to Eq. (|1^) involves 
only two approximations: (i) the frozen approximation for the configuration of the spectator 
system and (ii) the eikonal approximation for the trajectory of the struck particle. Explicit 
calculation of the relevant Mon(p,q;t) integrals within nonrelativistic many-body theory 
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i=2 



viri 



vr 



(13) 



appears to be feasible, at least for few-nucleon systems and infinite nuclear matter. However, 
to establish a clear connection with the PWIA picture, it is useful to devise approximations 
that permit ^(q, tu) of Eq. (|1^) to be expressed in terms of either the spectral function 
P(k, E) or the PWIA response >S'pwia(c[) f-i^)- 
The spectral function is defined as 

P(k, E) = Y^ \Mon{k)\'5{E + - E„) . (18) 

n 

To recover this quantity within the formula (15), the integrals Mq^ of Eq. ( |I7D are required 
to take the form 

Mon(p, q; t) = Mon(p - q) W(p, q; t) , (19) 

where the function U is to be independent of the state of the spectator system, labeled by 
the index n. An even more drastic simplification is achieved upon requiring that the time 
dependence of Mon(p,q;t) be factorizable according to 

Mon(p, q; t) = Mo„(p - q)f7q(t) , (20) 



i.e. upon assuming that the function U defined by Eq. ([19D does not depend upon p, which 
in turn corresponds to making the approximation p ~ q in W. 

Substitution of Eq. (pO|) into Eq. ([15|) allows one to rewrite the response as a convolution 
integral, 

S{q, uj)= J c?u;'5pwiA(q, oo')F^{uj - u') , (21) 

the folding function Fq{uj) being given by 

F^i^) = J ^^''''U^it) . (22) 

To obtain the function Uq{t) embodying all FSI effects {N.B. the PWIA can be regained by 
setting f/q(t) = 1, i.e., F{(jj) = 6{uj)), one makes the replacement 

Lrq(ri+vr-r,) ^ \ = ^q(r) , (23) 



J=2 



which amounts to averaging FSI with the ground-state configuration weights. In infinite 
nuclear matter at uniform density p, the average involved in Eq. (|23) takes the simple form 



V^{T)=pj£rgir)T^iv + VT), (24) 

where the radial distribution function g{r) measures the probability of finding two nucle- 
ons separated by a distance r = |r|. Keeping only the contributions associated with the 
imaginary part of the NN amplitude, which is known to be dominant at large incident 
momentum, we can finally write the eikonal propagator as 

f7q(t) = exp J drlmFq(r) . (25) 



The approach developed in I and employed in Refs. is based on the assumptions 

underlying Eqs. (|20D-(p5D. A different approximation scheme leading to the factorization of 
Mo„ can be obtained by inserting into Eq. (0) the identity 

/ dR'd\6{r, - r[) 5^(^-1) (i? -R') = J2 f ^R' ^*m{R')^n{R) , (26) 
J J 

where the sum includes a complete set of eigenstates of the A-particle Hamiltonian H. This 
procedure is not unique, because the 5-function insertion allows for different assignments of 



the arguments of the functions entering the right-hand side of Eq. (|T^ , leading to different 



but ultimately equivalent representations. Equating two such representations, we have 
^ / rfi2ciVie-^^'-i$;(^)^^(ri,i?) / rf^'dMe^'''^'"^"'''^^^o(r;,i?')n(ri,i?')f^pW,i?';t) 

= Y.[ dRd'ne-'''-''K{R)^N{R)Up{r,,R;t) [ dR'd'r[e'''<''^~<^^o{r[, R^A^u R') , 

(27) 

with k = p — q. Now imposing translation invariance, as appropriate for the case of infinite 
nuclear matter, and retaining only the term corresponding to |A^) = |0) in the above sums, 
Eq. ( pTl) becomes 

Mon(p, q; t) = Mo„(k) . (28) 

n(k) 
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Here n(k) is the nucleon momentum distribution, defined in terms of tlie spectral function 
as 

n(k) = J dEP(k,E), (29) 

while 

Uo{k, P;t) = J d'nr e''^-'-"'P(p, ru'; t) , (30) 

with rii/ = ri — ri/ and 

P(p,rn';t) = J dR%{r[,R)^o{ruR)Up{ruR;t) . (31) 

It can be readily seen that if Up{R, ri, t) = 1, i.e. if FSI are absent, the function V{p, rn'; t) 
reduces to p{ri]Ti')/A, where 

p^{r,-r[)=Aj dR%{r[,R)^o{ri,R) (32) 

is the one-body density matrix, whose Fourier transform is nCk). As a consequence, we have 
Mon(p,q;t) = Mo„(k) and the PWIA is recovered. In addition, the relation ( p8|) can be 
shown to hold as an equality at the lowest order of the cluster expansion. 
Substitution of Eq. ( pH]) into the definition of the response leads to 

5(q,a;) = Jduj'J -0^F^,^{uj' - io) J dE P{k, E)S{uj' - E - Ei^+^i) , (33) 

where the generalized folding function Fk,q(co') is defined as 

1 r dt 



Fi.A^) = ^J i^«"*^o(k, k + q; t) . (34) 



III. MANY-BODY CALCULATION OF THE GENERALIZED FOLDING 

FUNCTION 

Eqs. ( ppp and ( pi] ) show that calculation of the generalized folding function of Eq. ( |51[ ) 
requires a knowledge of the partially diagonal n-body density matrices 
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Pn{ri, r2, . . . , Tn] r[, r2, . . . , r„) = J^^^ J d^r^+i ■ ■ ■ d^rA%{r'i, i?)^o(ri, R) (35) 

of the target nucleus, for all n <A. 

The numerical calculation of p„ within an ab initio microscopic approach involves pro- 
hibitive difficulties, even for the case of infinite nuclear matter considered here. In view of 
this problem, we need to model p„ in terms of quantities that consistently incorporate the 
relevant physics and can still be reliably calculated. The results presented in this paper have 
been obtained using an approximation scheme (hereafter referred to as the hole approxima- 
tion), in which p„ is written in terms of the one-body density matrix of Eq. (^) and the 
half-diagonal two-body density matrix 

P2(ri, r2; r^, r^) = A{A - I) j d^r, . . . rfV^^*(r;, i?)*o(ri, R) . (36) 

The explicit formula for the resulting ra-body density matrix, 

HA / 1 " / 

Pn (i"!, r2, . . . , r„; r'^, r2, . . . , r„) = -f — — — ^ |j p2(ri, r^; r'^, r^) , (37) 

shows that the hole approximation represents the multinucleon spatial correlations involved 
in the definition of superposition of two-nucleon correlations. 

Among other expressions that can be constructed from the same building blocks, Eq. (|37D 
was chosen primarily because it fulfills some basic properties of the exact density matrices. 
In particular, p^^, which is obviously real, satisfies exactly the asymptotic factorization 
requirement 

lim p„(ri, r2, . . . , r^, r2, . . . , r„) = pp„_i(ri, r2, . . . , r„_i; r^, rg, . . . , r„-i) , (38) 

|rn|^oo 

while violating, although not severely, the sequential relation 

/dV„p„(r.r.....r„;r;.r......r„) 

= [A- (n- 1)] p„_i(ri, r2, . . . , r„_i; r^, r2, . . . , r„_i) . (39) 



Within the hole approximation, Eq. (|39D translates into 
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Tni 1"]^, r2, . . . , r„j 
= - l)Pn^i(ri, r2, • • • , rn-i; r'^, ra, 
= [A-(n-l)]p«^, + 0(n-2). 



(40) 
(41) 



In addition to p2, calculation of the generalized folding function in the hole approxima- 
tion calls for a knowledge of the imaginary part of the quantity Fq of Eq. (|1^), i.e. of the 
imaginary part of the A^A^ scattering amplitude /q. In this work, we have employed the 
simple parametrization originally proposed in Ref. namely 



Im/q(k) = ^CTTVArexp (-/3^|k|^) . 



(42) 



Numerical values of the total A^A^ cross section a^N and the slope parameter j3 resulting 
from fits to A^A^ scattering data are given in Refs. and . 

Using the hole approximation (^) and together with the parametrization (42) of ImFq, 
we can finally assemble the working expression 

1 



n(k) 



1 dz 
V Jo 2tt 



2Re- / ^ex.p{iuj-] / drw e'^""^^' Eq{z,rii'] 



(43) 



for the generalized folding function, with 



-Eq(2;,riiO 



exp 



<^NN 



(44) 



and 



2, 



, r dr dcosO P2irii> ,r,r') 
Piirw) Jo J-i 



X exp 



sin 6 

4/?2 



/2; + rcos6'\ fr cos 9' 

erf — erf 



(45) 



Iri — and r' = \r' 



213 ) V 2/9 

In the above equations, Piirni) and P2{,fiii ,r,r'), with r = |r| 
\t'i — ral, are the nuclear- matter one-body and half-diagonal two-body density matrices, 
respectively, v is the speed of the struck particle, the z axis is chosen along the direction 
of q, and 9 is the angle between r and q. It can be shown that the k-independent CGA 
folding function of Ref. can be recaptured from Eqs. (43)-(45) by replacing P2{ji, ^'i-, ^2) 



13 



with its fully diagonal part p2(ri7 ri, r2). The numerical calculation of Fk,q(co') has been 



performed by applying (i) the scheme proposed in Ref. |T8| to obtain n{k) and Pi(?"ii') 
and (ii) the formalism developed in Ref. |19[ to evaluate the nuclear-matter half-diagonal 
two-body density matrix. Both approaches use correlated many-body wave functions and 
FHNC integral equations to sum up selected cluster contributions to all orders in the relevant 
expansions. 

In particular, p2(ri, T2; r[, has been approximated by its leading dd part in the FHNC 



formalism according to thus we set 



P2(ri, ra; r^, ra) = pi{ru')9Qddir)gQdd{r') 



(46) 



where gQdd{f) consists of the dd nodal and non-nodal components of the FHNC expression for 
pi(ri; r'^). (The notation dd refers to the topological classification of the diagrams associated 
with the corresponding terms in the cluster expansion. See Ref. [0 for details.) 

Due to its weak dependence upon cos0 = f ■ rn/, the quantity gqddi^') is approximated 
by its angular average, according to 

1 



9Qdd{r[) gqddir') 



(i(cos0) gQdd {r^ + r^' - 2rr^^/ cos0)^/^ 



The resulting expression for the generalized folding function is 

4 dz 



^k,q(^) 



coslu;-) / r^ydriv jQ{krii,)pi{riii)E^{z,rii, 



n(k) Jo V \ vj Jo 
where jo{x) is the zeroth-order spherical Bessel function and 



E^{z, rir) = exp 



o'nn 
167r/?2 



Jq(2;,rii/) 



with 



00 

2, 



27rp r dr J^^dcosegQdd{r)gQdd{r') 
sin^ 9 



X exp 



4/?2 



jz + r cos 9\ ^(rcosO 
erf — err 



2P 



2/9 



(47) 



(48) 



(49) 



(50) 
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IV. NUMERICAL RESULTS 



We have carried out the calculation of the nonrelativistic response of a realistic model of 
nuclear matter, based on the Hamiltonian 

A 2 A A 

^ = E£^+ E E v^^^^ (51) 

i=l j>i=l k>j>i=l 

where Vij and Vijk are potentials describing two- and three-nucleon interactions, and on a 
variational ground-state wave function 

MR)=S[IlF^MR)- (52) 

In this trial form, 5 is a symmetrization operator acting on the product of two-nucleon 
correlation operators, Fij, and xo is the Slater determinant describing a noninteracting Fermi 
gas of nucleons with momenta k filling the Fermi sea, i.e., with |k| < kp = (37r^p/2)^/^. 
The operator Fij, which should reflect the correlation structure induced by the nuclear 



Hamiltonian, has been chosen as p2 



Fij = fc{rij) + fa{rij){(Ti ■ (Tj) + fr{rij){Ti ■ Tj) + faT{rij){cTi ■ (Tj){Ti ■ T j) 

+ ft{rij) Sij + ftr{rij) Sij {Ti -Tj). (53) 

Here Sij = 3(crj • rij){(Tj ■ rjj)/|r.yp — (cTj ■ cTj) is the usual tensor operator, while 
fc{r), fa{r), fr{r), far{r), ft{r), and ftrir) are correlation functions whose radial shapes are 
determined by minimizing the expectation value of the Hamiltonian (0) in the ground state 
described by Eq. (H) |2|. 



The PWIA response has been calculated using the nucleon spectral function of Ref. pOj 



obtained from a nuclear Hamiltonian including the Urbana f 14 A^A^ potential supplemented 
by the TNI model of the three-body interaction [^]. The same ingredients entering the 
calculation of the spectral function have been employed in the calculation of the density 
matrices needed to obtain the generalized folding function defined by Eqs. (^8D-(pO|). 

The dependence of the folding function Fk,q(c<j) on A; = |k| and u is illustrated in Figs. 1 
and 2 at |q| = 1.31 and 1.94 GeV/c, two values representative of the range covered by 
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the data analyzed in Refs. For comparison, we also show the /c-independent folding 

functions obtained using the approach developed in I. We note that in symmetric nuclear 
matter at its equilibrium density p = po = O.IQ fm~^, the Fermi momentum is kp = 1.33 
fm~^. It is apparent that in the region k < kp, the folding function Fk,q(i^) is very close to 
its CGA counterpart, whereas a strong ^-dependence is observed at k > kp. At |k| > kp 
the generalized folding function shrinks, and its tail begins to oscillate, implying that these 
k values correspond to weaker FSI. On the other hand, for larger momenta, well above the 
Fermi level, -Fk,q(i^) gets broader again. The results of Figs. |l] and ^ confirm the naive 
expectation that while the averaging procedure involved in the approach of Ref. is quite 
reasonable in the region of |k| < kp, where the nucleon momentum distribution is nearly 
constant, the momentum dependence of the folding function associated with fast nucleons, 
carrying momenta larger than kp, must be treated explicitly. 

Figs. ^ and | show the dependence of the nuclear matter y-scaling function 



evaluated for the respective choices |q| = 1.94 and 1.31 GeV/c of the momentum transfer. 
Using the y-scaling function rather than the response makes it easier to directly compare 
FSI effects at different values of |q|, as measured by the deviation from the PWIA results 
at fixed y. 

Comparison of the PWIA scaling functions (dash-dot lines) to the results obtained from 
the CGA (dashed lines) and the approach described in the previous sections (solid lines) 
clearly demonstrates that use of the k-dependent generalized folding function leads to a 
suppression of FSI effects in the region of large negative y, corresponding to the low-energy 
tail of the response [see Eq. (|55D]. For example, at y = 600 MeV/c the differences between 
the dashed and solid lines are ~ 70 % and ~ 20 % at |q| = 1.31 and 1.94 GeV/c, respectively. 



F{q,y) 




(54) 



on the scaling variable 




(55) 
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The fact that the suppression of FSI appears to be more pronounced at the lower values 
of |q| indicates that a broader CGA folding function is associated with a smaller effect of 
k-dependence. In fact, the larger value of the A^A^ cross section at |q| = 1.94 GeV/c, namely 
aNN — 43 mb compared to aNN — 35 mb at |q| = 1.31 GeV/c, makes the CGA folding 
function broader at the higher momentum transfer (see Figs.|I] and |^). 



V. SUMMARY AND CONCLUSIONS 

We have carried out a calculation of the nuclear matter response in which A^A^ correla- 
tions, which are known to play an important role in the low-energy tail of the response, have 
been taken into account both in the initial and final states. The effects of dynamical corre- 
lations in the initial state have been consistently incorporated into the PWIA calculation, 
based on a realistic spectral function obtained from an ab initio microscopic many-body 
approach [BO]. 



The same ingredients entering the calculations of the spectral function of Ref. |]20| have 
been employed in the calculation of corrections to the PWIA arising from FSI between the 
struck nucleon and the spectator system. Our many-body treatment of FSI is based on 
much the same scheme as applied in I and is therefore predicated on the eikonal and frozen 
approximations. However, the new treatment allows us to go beyond the simple convolution 
expression for the response and explicitly take into account the dependence on the initial 
momentum of the struck nucleon, which is averaged over in the CGA. 

Numerical results show that the momentum dependence of the generalized folding func- 
tion produces a sizable effect in the low-energy tail of the ?/-scaling function in the range 
of momentum transfer 1 < |q| < 2 GeV/c covered by the inclusive electron- nucleus scat- 
tering data analyzed in Refs. |@-|^. While repeating the analysis of Refs. [0-0] within the 
approach proposed here would certainly be of great interest, it must be pointed out that 
using a momentum-dependent generalized folding function would in no way help to improve 
the agreement between the CGA and the data. As shown in I, the discrepancy between 
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the CGA results and the measured cross sections increases as |q| increases, while the sup- 
pression of FSI due to the momentum dependence of the folding function appears to be 
larger at lower momentum transfer. A different mechanism leading to a quenching of FSI 
and exhibiting the opposite momentum-transfer dependence, such as the one associated with 
the color-transparency model employed in Refs. 0-^, still seems to be needed to reconcile 
theory and data. 
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FIG. 1. Dependence of the generalized folding function Fk,q(a;) of Eqs. ([4g|)- (|50|) on A; = |k| 
and uj, at q = |q| = 1.94 GeV/c. The diamonds show the results of the approach of Ref. in 
which the dependence of the folding function upon the initial nucleon momentum k is neglected. 
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FIG. 2. Same as in Fig|T|, but for momentum transfer ^7 = |q| = 1.31 GeV/c. 
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FIG. 3. Dependence of the nuclear-matter scaling function -F(q, y) defined by Eqs. (|5^)-(|55[) 
on the scaling variable y, at q = |q| = 1.94 GeV/c. The dash-dot line shows the PWIA result, while 
the solid and dashed lines correspond respectively to calculations carried out using the generalized 
folding function of Eqs. (^)-(|50|) and the |k|-independent CGA folding function of Ref. [^. 
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FIG. 4. Same as in FigJ^, but for momentum transfer g = |q| = 1.31 GeV/c. 
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